Above threshold ionization by few-cycle spatially inhomogeneous fields 
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We present theoretical studies of above threshold ionization (ATI) produced by spatially inhomo- 
geneous fields. This kind of field appears as a result of the illumination of plasmonic nanostructures 
and metal nanoparticles with a short laser pulse. We use the time-dependent Schrodinger equation 
(TDSE) in reduced dimensions to understand and characterize the ATI features in these fields. It is 
demonstrated that the inhomogeneity of the laser electric field plays an important role in the ATI 
process and it produces appreciable modifications to the energy-resolved photoelectron spectra. In 
fact, our numerical simulations reveal that high energy electrons can be generated. Specifically, 
using a linear approximation for the spatial dependence of the enhanced plasmonic field and with 
a near infrared laser with intensities in the mid- 10 14 W/cm 2 range, we show it is possible to drive 
electrons with energies in the near-keV regime. Furthermore, we study how the carrier envelope 
phase influences the emission of ATI photoelectrons for few-cycle pulses. Our quantum mechanical 
calculations are supported by their classical counterparts. 

PACS numbers: 42.65.Ky,78.67.Bf, 32.80.Rm 
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I. INTRODUCTION 

In the field of the interaction of laser fields with matter, 
above-threshold ionization (ATI) has been a particularly 
interesting subject in both experimental and theoretical 
physics. ATI, which was experimentally observed more 
than 30 years ago [l[, occurs when an atom or molecule 
absorbs more photons than the minimum number re- 
quired to ionize it, with the leftover energy being con- 
verted to the kinetic energy of the released electron. 

With recent advances in laser technology, it has be- 
come possible to generate few-cycle pulses, which find a 
wide range of applications in science, such as controlling 
chemical reactions and molecular motion 0, Q , and gen- 
erating high-order harmonics and even the creation of 
isolated extreme ultraviolet (XUV) pulses [3, H| . These 
allow even more control on an attosecond temporal scale. 

The electric field in a few-cycle pulse can be character- 
ized by its duration and by the so-called carrier-envelope 
phase (CEP). In comparison to a multicycle pulse, the 
electric field of few-cycle pulses is greatly affected by the 
CEP 0, 0] • The influence of CEP has been experimen- 
tally observed in high-harmonic generation (HHG) [8|, 
the emission direction of electrons from atoms M and in 
the yield of nonsequential double ionization [lOj. In or- 
der to have a better control of the system on an attosec- 
ond temporal scale it is, therefore, important to find reli- 
able schemes to measure the absolute phase of few-cycle 
pulses. 

Recently, the investigation of ATI generated by few- 
cycle driving laser pulses has attracted so much in- 
terest due to the sensitivity of the energy and angle- 
resolved photoelectron spectra to the absolute value of 
the CEP [ll|, [l2|. Consequently, this feature renders 
the ATI phenomenon a very valuable tool for laser pulse 
characterization. In order to characterize the CEP of 
a few-cycle laser pulse, the so-called backward- forward 



asymmetry of the ATI spectrum is measured and from 
the information collected the absolute CEP can be ob- 
tained [l3|| . Furthermore, nothing but the high energy 
region of the photoelectron spectra appears to be sen- 
sitive to the absolute CEP and consequently electrons 
with kinetic energy are needed in order to characterize 
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New experiments have demonstrated that the har- 
monic cutoff and electron spectra of ATI could be ex- 
tended further by using plasmon field enhancement [Tf| 
This field appears when a metal nanostructure or 
nanoparticle is illuminated by a short laser pulse and it is 
not spatially homogeneous, due to the strong confinement 
of the plasmonics spots and the distortion of the electric 
field by the surface plasmons induced in the nanosys- 
tem. One should note, however, that a recent controversy 
about the outcome of the experiments of Ref. [16| has 
arisen [18-20]. Consequently, alternative systems to the 
metal bow-tie shaped nanostructures have appeared [2l[ . 
A related process employing solid state targets instead of 
atoms and molecules in gas phase is the so called Above 
Threshold Photoemission (ATP). This laser driven phe- 
nomenon has received special attention recently due to its 
novelty and considering new physics could be involved. 
In ATP electrons are emitted from metallic surfaces or 
metal nanotips and they present distinct characteristics, 
namely higher energies, far beyond the usual cutoff for 
noble gases and consequently the possibility to reach sim- 
ilar electron energies with smaller laser intensities (see 
e.g. |22l-[27|). Furthermore, the photoelectrons emitted 
from these nanosources are sensitive to the CEP and 
consequently it plays an important role in the an gle and 
energy resolved photoelectron spectra [13, HH l28l l29| . 

Despite new developments, all numerical and semiclas- 
sical approaches to model the ATI phenomenon are based 
on the assumption that the external field is spatially ho- 
mogeneous in the region where the electron dynamics 
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take place [30ll3lj|. For an inhomogeneous field, however, 
important changes will occur to the features of strong 
field phenomena [H, [I?} since the laser-driven electric 
field, and consequently the force applied to the electron, 
will also depend on position. Up to now, there have been 
very few studies to investigate the strong field phenom- 
ena in such kind of fields |32i-l35j. 

From a theoretical viewpoint, the ATI process can be 
tackled using different approaches (for a summary see 
e.g. [l4l l36l440j and references therein). In this article, we 
concentrate our effort in extending one of the most and 
widely used approaches: the numerical solution of time- 
dependent Schrodinger Equation (TDSE) in reduced di- 
mensions. We have developed our numerical tool in such 
a way to allow the treatment of a very general set of non- 
homogeneous fields. Furthermore, based on our model, 
we examine the influence of the CEP on photoelectron 
spectra of ATI. The kinetic energy for the rescattered 
electron is classically calculated and compared to our 
quantum mechanical approach. 

This article is organized as follows. In Sec. II, we 
present our theoretical approach to model ATI produced 
by nonhomogeneous fields. Subsequently, in Sec. Ill, we 
employ this method to compute the ATI energy-resolved 
photoelectron spectra using few-cycle laser pulses for 
both homogeneous and inhomogeneous fields. In ad- 
dition, we perform classical simulations to support our 
quantum mechanical method. Finally, in Sec. IV, we 
conclude with a short summary and outlook. 



II. THEORETICAL APPROACH 

In order to calculate the energy resolved photoelec- 
tron spectra, we use the one-dimensional time-dependent 
Schrodinger equation (1D-TDSE) 



1 dt 



1 d 2 
2d^ 



+ Vatom(x) + Vl aser 



(1) 



where Vi aser (x,t) represents the laser-atom interaction. 
For the atomic potential, we use the quasi-Coulomb or 
soft core potential 



V a torn{%) — 



(2) 



which was introduced in [4JJ and has been widely used 
in the study of laser-matter processes in atoms. The 
parameter a in Eq. (|2j) allows us to match the ionization 
potential of the atom under consideration. We consider 
the field to be linearly polarized along the x-axis and 
modify the interaction term Vi aser (x,t) in order to treat 
spatially nonhomogeneous fields, although maintaining 
the dipole character. Consequently we write 



where E(x,t) is the laser electric field defined as 

E(x, t) = E f(t) (1 + eh(x)) sin(ut + 0). (4) 

In Eq. (J2J), Eo, uo and are the peak amplitude, the fre- 
quency of the laser pulse and the CEP, respectively. We 
refer to sin(cos)-like laser pulses where = (0 = tt/2). 
The pulse envelope is given by f(t) and e is a small pa- 
rameter that characterizes the inhomogeneity strength. 
The function h(x) represents the functional form of the 
nonhomogeneous field and, in principle, could take any 
form and be supported by the numerical algorithm |35j . 
In this work, however, we concentrate our efforts on the 
simplest form for h(x), i.e. the linear term: h(x) = x. 
This choice is motivated by previous investigations in 
high-order harmonic generation 

In the linear model we are using in this work, the units 
of e are inverse length (see also |32l434j). We have writ- 
ten Viaser hi Eq. (|3]) in such a way to emphasize the 
fact we are working within the dipole approximation and 
any deviation of it is considered small, i.e. higher electric 
multipole terms and magnetic effects are neglected (43| . 
To model short laser pulses, we use a sin-squared enve- 
lope f(t) of the form 



/(*) 



sin 



cot 
2n n 



(5) 



Vl aser (X,t) = -E(x,t)x 



(3) 



where n p is the total number of optical cycles. The total 
duration of the laser pulse will then be T p = n p r where 
r = 2tt/u is the laser period. 

We assume the target atom is in the ground state (Is) 
before we turn on the laser (t = — oo). This state can be 
found by solving an eigenvector and eigenvalue problem 
once the spatial coordinate x has been discretized. We 
chose a 2 = 1.412 to match the atomic ionization potential 
of our target, which is an hydrogen atom (I p = —0.5 
a.u.). Eq.(pQ) is solved numerically by using the Crank- 
Nicolson scheme with an adequate spatial grid [30]. We 
employ boundary reflections mask functions [44| in order 
to avoid spurious contributions. 

For calculating the energy-resolved photoelectron spec- 
tra P(E) we use the window function technique devel- 
oped by Schafer [45l. 146] . This tool has been widely used, 
both to calculate angle-resolved and energy-resolved pho- 
toelectron spectra [47j and it represents a step forward 
with respect to the usual projection methods. 



III. RESULTS 

In this section, we will determine the energy-resolved 
photoelectron spectra P(E) using Eq. (pQ), in order to 
investigate the role of the inhomogeneities of the field. 
Furthermore, we demonstrate how the CEP will effect 
the the energy-resolved photoelectron spectra of ATI. We 
employ a four-cycle (total duration 10 fs) sin-squared 
laser pulse with an intensity I = 3 x 10 14 W/cm 2 and 
wavelength A = 800 nm. 
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We chose three different values for the parameter that 
characterizes the inhomogeneity strength, namely e — 
(homogeneous case), 0.003 and 0.005. Figures 1 and 2 
show the cases with <j) = (a sin-like laser pulse) and 
(j) = 7r/2 (a cos- like laser pulse), respectively. Panels 
(a) of both Figures represent the homogeneous case, i.e. 
e = 0, and panels (b) and (c) show the nonhomogeneous 
case with e = 0.003 and e = 0.005, respectively. 
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FIG. 1. Energy-resolved photoelectron spectra P(E) calcu- 
lated using the 1D-TDSE for a model atom with I p = —0.5. 
The laser parameters are I = 3 x 10 14 W/cm 2 and A = 800 
nm. We have used a sin-squared shaped pulse with a total 
duration of 4 cycles (10 fs) and (j) = (a sin- like pulse). The 
arrows indicate the 2U P and 10U P cutoffs predicted by the 
classical model 14]. Panel (a) £ = (homogeneous case), (b) 
e = 0.003 and (c) e = 0.005. 



For the homogeneous case, the spectra exhibits the 
usual distinct behavior, namely the 2U P cutoff (~ 36 
eV for our case) and the 10U P cutoff (« 180 eV), where 
Up = Eq/Auj 2 is the ponderomotive potential. The for- 
mer cutoff corresponds to those electrons that, once ion- 
ized, never return to the atomic core, while the latter 
one corresponds to the electrons that, once ionized, re- 
turn to the core and elastically rescatter. It is well es- 
tablished using classical arguments that the maximum 
kinetic energies of the direct and the rescattered elec- 
trons are E^^ = 2U P and E^^ = 10Z7 P , respectively. 
In a quantum mechanical approach, however, it is pos- 
sible to find electrons with energies beyond the 10/7 P , 
although their yield of them drops several orders of mag- 
nitude [14| . Experimentally, both mechanisms contribute 



to the energy-resolved photoelectron spectra and conse- 
quently the theoretical approach to tackle the problem 
should to include them. In that sense the TDSE, which 
can be considered as an exact approach to the problem, 
is able to predict the P(E) in the whole range of electron 
energies. In addition, the most energetic electrons, i.e. 
those with Ek ^> 2U P , are used to characterize the CEP 
of few-cycle pulses. As a result, a correct description of 
the rescattering mechanism is needed. 
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FIG. 2. Idem Fig. 1 but <j) = tt/2 (a cos- like pulse). 

For the inhomogeneous case, the cutoff positions of the 
direct and the rescattered electrons are extended towards 
larger energies. For the rescattered electrons, this ex- 
tension is very prominent. In fact, for e = 0.003 and 
e = 0.005, it reaches « 260 eV and w 420 eV (panels 
b and c of Fig. 1, respectively). Furthermore, it ap- 
pears that the high energy region of P(E), for instance, 
the region between 200 — 400 eV for e = 0.005 (see pan- 
els (c) of Figs. 1 and 2), is strongly sensitive to the 
CEP. This feature indicates that the high energy region 
of the photoelectron spectra could resemble a new and 
better CEP characterization tool. It should be, how- 
ever, complemented by other well known and established 
CEP characterization tools, as, for instance, the forward- 
backward asymmetry (see [14(). Furthermore, the uti- 
lization of nonhomogeneous fields would open the avenue 
for the production of high energy electrons, reaching the 
keV regime, if a reliable control of the spatial and tem- 
poral shape of the laser electric field is attained. 

We now concentrate our efforts in order to explain 



4 



8 



^ 200 




(a) 



i 



(b) 



(c) 




Direct electrons 
Rescattered electrons 



(d) 



(e) 



(0 




3 4 1 2 3 4 

Time (in optical cycles) 

FIG. 3. (Color online). Numerical solutions of the Newton equation (Eq. ([6])) plotted in terms of the direct and rescattered 
electron kinetic energy. The laser parameters are the same as in Fig. 1. Panels (a), (b) and (c) correspond to the case of sin- like 
pulses (0 = 0) and for e = (homogeneous case), e = 0.003 and e = 0.005, respectively. Panels (d), (e) and (f) correspond to 
the case of cos- like pulses (<fi = 7r/2) and for e = (homogeneous case), e — 0.003 and e — 0.005, respectively. 



the extension of the energy-resolved photoelectron spec- 
tra using classical arguments. From the simple-man's 
model [48| we can describe the physical origin of the ATI 
process as follows: an atomic electron at a position x = 0, 
is released or born at a given time, that we call ioniza- 
tion time U, with zero velocity, i.e. x(U) = 0. This elec- 
tron now moves only under the influence of the oscillat- 
ing laser electric field (the residual Coulomb interaction 
is neglected in this model) and will reach the detector 
either directly or through the rescattering process. By 
using the classical equation of motion, it is possible to 
calculate the maximum energy of the electron for both 
direct and rescattered processes. The Newton equation 
of motion for the electron in the laser field can be written 
as j3]): 



electron released or born at time ti is 



x(t) 



-VxVla 

E(x,t) 
E(t)(l- 



er £) 

- [V x E(x,t)]x 
2sx(t)), 



(6) 



where we have collected the time dependent part of the 
electric field in E(i), i.e. E(t) = E f(t) sm(ut + aV) and 
we have specialized to the case h(x) = x. In the limit 
where e = in Eq. (j6]), we recover the homogeneous 
case. For the direct ionization, the kinetic energy of an 



E d 



[x(ti) - x(t f )Y 



(7) 



where tf is the end time of the laser pulse. For the rescat- 
tered ionization, in which the electron returns to the core 
at a time t r and reverses its direction, the kinetic energy 
of the electron yields 



E r 



[x(U) + x(t f )-2x(t r )Y 



(8) 



E, 



For homogeneous fields, Eqs. dH and dU become as 

= [A { n)-A {tf) f ^ = [AiUHAit^AiU)]^ with 

A(t) being the laser vector potential A(t) = — J 1 E(t')dt' . 
For the case with e = 0, it can be shown that the max- 
imum value for Ed is 2U P while for E r it is 10U P [H|. 
These two values appear as cutoffs in the energy resolved 
photoelectron spectrum as can be observed in panels (a) 
of Figs. 1 and 2 (see the respective arrows). 

In Fig. 3, we present the numerical solutions of Eq. 
([6]) , which is plotted in terms of the kinetic energy of the 
direct and rescattered electrons. We employ the same 
laser parameters as in Figs. 1 and 2. Panels (a), (b) and 
(c) correspond to the case of <p = (sin-like pulses) and 
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for e = (homogeneous case), £ = 0.003 and £ = 0.005, 
respectively. Meanwhile, panels (d), (e) and (f) corre- 
spond to the case of <p = tt/2 (cos-like pulses) and for 
e = (homogeneous case), e — 0.003 and e — 0.005, re- 
spectively. From the panels (b), (c), (e) and (f) we can 
observe the strong modifications that the nonhomoge- 
neous character of the laser electric field produces in the 
electron kinetic energy. These are related to the changes 
in the electron trajectories (for details see e.g. (33l-[35|). 
In short, the electron trajectories are modified in such a 
way that now the electron ionizes at an earlier time and 
recombines later, and in this way it spends more time 
in the continuum acquiring energy from the laser electric 
field. Consequently, higher values of the kinetic energy 
are attained. A similar behavior with the photoelectrons 
was observed recently in ATP using metal nanotips. Ac- 
cording to the model presented in Ref. [27] the localized 
fields modify the electron motion in such a way to allow 
sub-cycle dynamics. In our studies, however, we consider 
both direct and rescattered electrons (in Ref. [27| only 
direct electrons are modeled) and the characterization 
of the dynamics of the photoelectrons is more complex. 
Nevertheless, the higher kinetic energy of the rescattered 
electrons is a clear consequence of the strong modifica- 
tions of the laser electric field in the region where the 
electron dynamics takes place, as in the above mentioned 
case of ATP. 



IV. CONCLUSIONS AND OUTLOOK 

We have extended previous studies of high-order har- 
monic generation produced by nonhomogeneous fields 



to above threshold ionization (ATI). An example is the 
field generated in a vicinity of a metal nanostructure or 
nanoparticle when it is irradiated by a short laser pulse. 
We have modified the time dependent Schrodinger equa- 
tion to model the ATI phenomenon driven by nonhomo- 
geneous fields. We predict an extension in the cutoff po- 
sition and an increase of the yield of the energy-resolved 
photoelectron spectra in certain regions. These features 
are reasonable well reproduced by classical simulations. 
Our predictions would pave the way to the production of 
high energy photoelectrons, reaching the keV regime, us- 
ing plasmon enhanced fields. Application of our model to 
a broader range of laser parameters, including an exhaus- 
tive study of CEP effects, and a systematic survey over 
different atomic species using a full dimensional scheme 
will be subject of further investigations. 
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